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Abstract. 

Dub rouLI ne, oTADET' analyzes a Structure co iindsit 
it is stable and determinate. The structure must first be 
idealized as a rigid body supported by a system of links. 
The program checks on the arrangement and number of the 
Support links. Flags are set so that this subroutine may 
be used to select whether to use simple statics or strength 
of materials methods in subsequent analysis of the problem. 
Testing of the subroutine is discussed. 

Two input-output programs are also described. One 
calculates the reactions in the equivalent support linkages 
from the resultant load; the other assumes that the pro- 
perties of the links are known and calculates the allowable 


oad. 
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KPROB 
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A matrix of coefficients 


The coefficients of the equation of a 
plane 


A vector 

Rectangular components of A 
A column matrix 

A vector 

Rectangular components of B 


Vectors made up of the direction cosines 
of the support links, link cosine vectors 


Indicates the ith, jth or kth of B,, B 


? 
euc. 2 


A Fortran floating point, numerical 
format 

A vector force 

The ith vector force 

Rectangular components of F 
Rectangular components of F; 

A Fortran integer, numerical format 


Unit orthogonal base vectors in the 
X, y, Z2 directions respectively 


An array used to indicate the kind of 
problem involved 


Direction cosines of a line 
Vector moment 
Rectangular components of M 


An array used to indicate the reason 
for the lack of a solution 


V1 


ae Po» etc. 


PLA 


Position vectors giving the locations 

of the points of application of the 
links with respect to the first point 

of applicat@on 

An array used to store the coefficients 
of the equations of the planes determin- 
ing the lines of action of the links 


The position vector of the point of 
application of a force 


Rectangular components of r 
Poynting's vector 

A weight 

A column matrix 

Rectangular coordinates 
Coordinates of a given point 


The scalar result of forming a dot 
prodquct 


The vector result of forming a cross 
product 


The angle included between two vectors 
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ieee antroduction. 

Digital computers have been extensively applied to 
the solution of mechanics problems. However, these ap-~ 
plications have been limited in scope. The computer has 
been used as a "Super slide rule" to solve a particular 
type of problem which was too time-consuming or too dif- 
ficult previously. It is the intention of Professor Gilles 
Cantin to attack the problem of the design of structures 
by a digital computer in a general rather than a specific 
manner. 

One of the first routines needed for the Cantin Pro- 
ject was a routine to reject unstable structures and, 
perhaps, make some selection as to the methods of analysis 
to be used in subsequent tests. The subroutine "STADET" 
analyzes an idealized structure to find if it is stable 
and determinate. It sets flags which can be used to print 
out a description of the kind of problem involved. When 
the subroutine finds the first condition which would pre- 
clude the solution of the problem by simple statics, a 
flag is set and control returns to the input-output program. 
If the "error" is caused by a redundant support system, 
other methods of analysis may be available whereas there 
is no solution for an "unstable" problem. 

The definition of stability which is considered here 
is more restrictive than the question of stable equili- 
brium usually encountered in basic statics courses. Here 


is meant, not the criterion of energy change with small 





perturbations, but the geometrical stability of infinitely 


strong members. As expressed by Kinney, 1 


A stable structure will support any conceivable 
system of applied loads, resisting these loads 
elastically and immediately upon their applica- 
tion the strength of all members and the capacity 
of all supports being considered infinite. 


Consider the two links supporting a weight shown in 


Fee =). 





Fig. 1. Support links before and after displacement. 


In Fig. la, with the links horizontal, there is no 
vertical component of the tension in the bars to support W. 
Even with infinitely strong foundations and members, the 
geometrical arrangement precludes a vertical component une- 
til some downward displacement of W has occurred (Fig. 1b). 
The methods of strength of materials are required to find 
mae equilibrium position of the structure. Since W is not 
resisted "immediately upon its application," this arrange- 
ment fails the test for stability. 

Both a structure and its support system must be 
stable. The problem which is examined here is limited 


to that of the external stability of a structure. The 


aie Sterling Kinney, Indeterminate Structural Analysis, 
Pp. 20, oD 7 





structure itself 1s considered to be a rigid body so that 
only the system of supporting links is of concern. No 
loads are imposed during the stability analysis since 
equilibrium is not involved. 

However, in some cases forces such as gravity must be 
considered in order to have a structure which meets this 
Restricted definition of stability. This is a familiar 
assumption in basic statics where, for example, the gravity 
force provides a reaction against upward motion of a roller 
mest. it is up to the user of this program to decide 
whether a "gravitational" restraint is sufficient or whether 
a positive restraint is needed. 

The user must first idealize his structure so that the 
subroutine can operate on it. Appendix A, which shows the 
linkage diagrams of some common supports, may be helpful in 
this regard. However, it is at this point that engineering 
judgment is required. Needless to say, if the proper equi- 
valent support system is not employed, no reliance can be 
placed on any answers which are produced. 

Fig. 2 shows the decisions which are made by the sub- 
routine. The titles within the blocks indicate the basic 
sections of the subroutine. These titles are roughly indica- 
tive of the purpose of the sections though more information 
is usually produced than the title would lead one to belisve. 

For input the subroutine requires the direction cosines 
of the support linkages, the coordinates of the points of 


application of the links and the number of links. Appendix H 





pigs 2 oimplified Flow Chart 


COPLANARIL 
Reunes of application of links. coincident 
No Yes Unstable 
Points of application on a line? <3 Unstable 
No Yes 3 links?|]>3 Redundant 
Yes 
Points of application coplanar? 
No Yes 
,Links in same plane? 
No Yes 
<3 Unstable 
3 links? 1>3 Redundant 
Yes 


<6. Unstable 
6 links? |1>6. Redundant 





‘Links parallel? Links parallel? 
No Yes. Unstable No Yes 
Unstable 
Links in parallel planes? 
No Yes Unstable Links coplanar: 
Yes {No 


PARALLELZ Unstable 


Of Of links. parallel 
No Yes. Unstable 


CONCURRENT 


Lines of action of links intersect? 
No Yes Unstable 


y¥Structure is satisfactor 








shows how the input must be arranged to be read in by the 
input-output programs of Appendix ©. The output consists 
of the two flag arrays, KPROB and NOSOL. 

The subroutine first considers the points of appli- 
cation of the supporting links. If these points are along 
a line, the question of whether it is a two- or three- 
dimensional problem depends on the orientation of the links. 
If the points are in a plane, the links must all be in that 
Same plane for a two-dimensional problem to be involved. If 
the points are randomly located in space, only a three- 
dimensional problem can be involved. The two-dimensional 
problem requires three linkages for stability; the three- 
dimensional problem, six. If there are fewer linkages, an 
"unstable" error print-out is made; if more, an "indeter- 
minate" one is made. 

When it has been determined that the correct number of 
support linkages is involved, the orientation of the links 
is considered. The following restrictions apply to the 
three-dimensional case: 

All the links may not be parallel. 

They may not lie in parallel planes. 

The axes of the links may not intersect one 

straight line since limited rotation would 

be possible about that line. 

4. The axes of the links may not intersect in a 
point or limited rotation is possible about 
that point. 

5. Four of the six links may not be parallel or 
intersect in any oné point for, in these cases, 
a straight line can be found through the point 


(or infinity) which intersects the axes of the 
two remaining links and restriction 3. is violated. 


WNFH 


In the two-dimensional case only restrictions 1. and 


Ee epoyosllniae 








Where the links are all parallel to a given line, 
limited motion is possible in a direction perpendicular 
to this line. 


In Fig. 3, where the pairs of links are in parallel 


Fig. 3. Arrangement with links in parallel planes. 


planes, some motion is possible in a direction perpen-= 
dicular to these planes before a restoring component is 


developed. In Fig. 4 some rotation about the intersection 


e—— 5} - ~ = 
| 
PS 
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et 


[hale Oe Arrangement with intersecting lines of action 
of links. 


Lp 


of the lines of action of the links is necessary before 
a restoring couple is developed. 
The following section and the appendices describe 


the "STADET" subroutine, its auxiliary subroutines and 
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the Input-Output Programs which were used to test "STADET". 
Fortran language was used for all the programs. The testing 


was done on a Control Data Corporation 1604 digital computer. 









2. Method. 

Subtitles indicate the sections of the flow charts 
(Figs. 5 through 9) and the "STADET" list (Appendix J) 
under discussion. Circled numbers on the flow charts are 
used for two purposes. The first is to facilitate the 
comparison with the list of Appendix J; the second is to 
indicate a transfer to another section of the program 
when the use of a line to show the path would be cumber- 
some. 

The subroutine only operates up to the detection of 
the first error. The appropriate flag is then set anda 
return is made to the Input-Output Program. 

a. COPLANARL 

This section of the subroutine generates the 

position vectors which give the locations of the points 

of application of the links. It checks to see whether 
these points are coincident, lie on a line, are ina plane 
or are randomly located in space. In the case where the 
points of application are in a plane, this section of the 
program checks to see if the links are in the same plane 
and so determines whether a planar problem is involved. 

After deciding whether a two- or three-dimensional 
problem is involved, this section counts the number of links 
to find out if the number is sufficient for stability. 

The computer operates with a binary approximation to 
a decimal number when it is used in the floating point 


mode. Hach calculation is approximate and round-off errors 





Fig. 5 Flow Chart of COPLANAR1 


START 
DEFINE ZERO 


GENERATE P's 


Note: At least one P must be non-zero 
or STADET will not operate. 









<3 (42) NOSOL (3 )=1}-JRETURN 
2 (42) --{KPROB(1)=1} ~[KPROB(L)=1 @ 


>3 (43) NOSOL(4)=1 RETURN 


ALL 
CROSSED? 
Note: 
Pi//P, 
YES 


€ 3 


KPROB(2)=1 NBAR =? (42) 


KPROB(4)=1 


Note: j # i, j#k 





<6r+(6)-~{u0S0L (3 )=1}- —~{RETUAH 


$-(84)--{G0 10 PARALLELI}~(76) 


YES > 6 Les} fNOSOL(4 )=1} {RETURN 
NO ALL 
DOTTED? 
YES 


(441) <2 +3 (3) 





3 
KPROB(4)=1}—-4Go TO PLANARI (74) 





can accumulate. For this reason the “zero” which must. be 
used in the program is not the mathematical zero but an 
allowable tolerance. The first portion of this section 
examines the coordinates of the points of application and 
determines the value of "zero" to be used throughout 
"STADET". Appendix L discusses the method. 

Vectors are generated between the first point of 
application and the other points. These vectors are ex-~ 
amined to find the first non-zero one. If none of them 
is non-zero, there is only one point of support, the struc- 
ture is unstable, and a return is made to the Input-Output 
Program. 

If a non-zero vector is found, the subroutine deter- 
mines whether it is the last one or not. If it is not 
the last, the other vectors are examined to find whether 
any others are non-zero. If the first non-zero vector is 
the last of the array of position vectors, or if all after 
the first are zero, the points of application of the links 
lie on one line. 

Having determined that the points of application are 
on a line, the subroutine considers the number of links. 

A three-dimensional linkage system is unstable because the 
links intersect a line and only a planar problem has a 
unique solution. Therefore, if there are more than three 
links, the system is redundant. If there are less than 
three links, the system is unstable. If the number of 


links is not satisfactory, the appropriate flag is set. 


LO 





Where a second non-zero vector, P5, is found, its 
cross product with the first, Pa; is formed. The cross 
product vector, 4,, 1s examined to find out if it 1s zero. 
If it is zero, P, and Po are parallel. Since P, and Po 
have the same initial point, this mades them coincident 
and indicates that the first three points of application 
are on a line. The subroutine then checks to see if all 
of the position vectors have been examined. If they have 
not all been examined, the search for another non-zero 
vector continues. If they all have been examined and all 
of the points of application are found to lie on one line, 
an investigation of the number of supporting links is made, 
as above. 

If a cross product vector, as found above, is not 
zero, a check is made of the number of supporting links. 
If there are only three links, there are only two position 
vectors which (since they are not parallel) determine a 
plane. If there are more than three links, there are more 
than the two vectors that were crossed to form are To 
check whether the points of application are ina plane, 


it is necessary to dot these other vectors with X,. The 


douse product, X a scalar, is compared with zero. If all 


1? 
the dot products are zero, the points of application lie 
in one plane. 

With the points of application of the links in a plane 


a stable system is possible with either a two- or three- 


dimensional linkage system. If all the links lie in the 


aEsE 





plane that contains the points of application of the links, 
the problem is one in two dimensions. To check for a two- 
dimensional problem, the vectors made up of the direction 


cosines of the links are dotted with xX. and the dot products 


i 
are compared with zero. Since the link coSine vectors have 
one point in the point-of-application plane, a zero dot 
product indicates that the vector lies entirely in that 
plane. If a planar problem is what is involved, i.e., all 
the dot products are zero, the proper number of links is 
checked for, as before. 

When the points of application are not in a plane or 
when the links are not in the point-of-application plane, 
a three-dimensional problem is involved and six links are 
required for stability. Suitable error flags are set if 
there are not six supports. 

At this point it has been determined whether a two- 
or three-dimensional problem is involved and that the 
correct number of supporting links is involved for the 
type of problem. Next, instability due to incorrect 
orientation of the links is investigated. 

b. PARALLEL 

This section checks to see if the supporting 
links are parallel or lie in parallel planes. The first 
link cosine vector is crossed with each of the succeeding 
ones and the resulting cross product vectors are compared 
with zero. If all of the cross products are zero, all of 


the links are parallel. 


IZ 








Fig. 6 Flow Chart for PAKALLEL1 





COUNT 1 






ALL 
DOTTED ? 


YES 
a _TES yf 506) NOSOL(6)=1 }®4 RETURN 
NO 


| Go LO PARALLELa}+( 38 ) 
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NOSOL(2)=1}- >{k ET URN 








The second portion of this section is concerned with 
discovering whether a situation similar to the one shown 
in Fig. 3 exists or not. In the above test a non-zero 
cross product will result if two ot the links are not 
wemallel. (If the situation is like that of Fig. 3, the 
Saoss product vector is in the direction of motion of the 
rigid body.) If all of the links are perpendicular to 
this cross product vector, the case is indeed similar to 
that of Fig. 3. This section checks whether the links 
are perpendicular by dotting each of the link cosine 
vectors with the cross product vectors and comparing the 
mesultant scalar with zero. if all of the dot products 
are zero, all of the link cosine vectors are in parallel 
planes. If all are in parallel planes, the appropriate 
flag is set; if not, the program goes on to the next 
section. 

ce PARALLEL2 

Even though all of the links are not parallel, 
the supporting linkage system is unstable, by criterion 
3., if four of the six links are parallel. Therefore, 

a separate check is made for this possibility. The cross 
product test is used along with a counting system. If the 
linkage system passes this test, the remaining check for a 
three-dimensional system is done by CONCURRENTI. 

d. PLANARL 

This section of the subroutine checks the three-~ 


link cases to see if the links are in a plane and if the 


14 








Fig. 7 Flow Chart for PARALLEL2 






1 COUNT 1 


ae 


ALL 
CROSSED? 


YES 


COUNT 26 >—TES 


NO 
GO TO CONCURRENTL 


iD 





Fig. 8 Flow Chart of PLANARL 





| YES 
NO (05) 4NOSOL(2)=i}--—>RETURN 
| Note: All three B's parallel 


YES 





GO TO CONCURRENT1| 
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links are parallel. In Fig. 2 there are two paths from 
COPLANAR1 to PLANARL. The first of these three-link cases 
has the points of application on a line. In the other 
case the points of application are in a plane and it is 
known that the links are in the same plane. For this 
second case there is some duplication of testing, but 

the extra operating time is negligible and there was some 
convenience in programming this way. 

This section crosses the first link cosine vector 
misc tne second. if the cross product is zero, the first 
1S crossed with the third link cosine vector. If this 
cross product is also zero, the three links are parallel 
and the structure is unstable. If the second cross product 
ms mot zero, it is dotted with the second link cosine 
fweewor, tf the dot product 1s zero, all three link. cosine 
vectors are in the same plane, a stable planar problem is 
possible and the program goes to CONCURRENTL. If the dot 
product is not zero, the three links are not coplanar and 
the system is unstable. 

When the first cross product is not zero, it is dotted 
with the third link coSine vector and the dot product is 
compared with zero. If the dot product is not zero, the 
System is unstable, as above. If the dot product is zero, 
the program goes to CONCURRENTI. 

e. CONCURRENTL 

The remaining check is to find if the lines of 


action of the links intersect at a point. Jn order to 


ly/ 





Fig. 9 Flow Chart of CONCURRENTL 
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NO | YES 
SOLVE [RETURN] 
TRIAD 7 
SUBSET 
NO 
ALL 
NO SUBSETS\. YES 
CONSISTENT? \—s¢ EXAMINED? > , 
Note: The two lines 
are parallel. 
| YES 
~~ 


PLANAR ONS 


6 E « \ 188 | 
~~ ee hes -[NOSOL{1)=1} +RETURN] 


{ NO NO 
202 


No/8s OR MORE YES , 
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be able to operate mathematically, it is nécessary to 
Convert the link cosine vectors to lines in Space. The 
line form which was chosen was that of the intersection 
of two planes. (See Appendix G for discussion of sub- 
routine PLANE.) 

All of the planes are determined and the coefficients 
are stored in the array PLA. The equations in PLA are 
examined four at a time but with an overlap, i.e., equations 
one through four are examined first, then equations three 
through six; and, if a six-link problem which gives 12 
equations is involved, the process continues until all 
12 have been examined. 

The three=-link system must have at least one inter-~ 
section or the links would all be parallel, a possibility 
which has previously been eliminated. In a six=link 
problem at least four of the lines of action of the links 
must intersect to make the linkage system unstable. There-~ 
fore, the system of examining pairs of lines whose equations 
are adjacent in the array PLA must find the invalidating 
intersection and a more complicated checking system is not 
eequLred . 

The coefficients of four equations are treated as 
though they were the coefficients of four simultaneous 
equations in four unknowns. They undergo the matrix 
Singularity test of REAC] (See Appendix F). If the de- 
terminant of the matrix is not zero, the set of equations 


is consistent and the first two lines do not intersect. 
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If the next overlapping set of four is also consistent, 
three lines do not intersect and the link support system 
is satisfactory. This would complete the checks by "STADET". 

On the other hand, if the determinant is zero, the set 
of equations is dependent or intonsistent and the two lines 
intersect or are parallel. If the lines intersect, the 
point of intersection can be found by solving three of the 
equations simultaneously. If the lines are parallel, the 
four combinations of these four equations, taken three at 
a time, are all inconsistent. No account is taken of this 
parallelism since it is already known that less than four 
links are parallel. The next overlapping set of four equa- 
tions is then considered. 

When an intersection is found, the coordinates are 
substituted into the equations of all the planes. The 
number of equations which are satisfied is counted. In 
the case of a planar problem, if six are satisfied, all 
the lines of action intersect and the linkage system is 
unstable. In the three-dimensional case, if eight or 
more are satisfied, the linkage system is unstable. Where 
instability exists, the appropriate flag is set. 

This completes the checks made by "STADET". The 
rigid body is supported by a system of links whose re-~ 
actions can be determined by simple statics or the appro- 
priate flag has been set in NOSOL. 

The one remaining possible error, that of not having 


the load in the correct plane for a planar problem, is not 


a0 








strictly of concern in deciding whether the system of 
links is stable and determinate. The check for this 


error is discussed in Appendix C, Input-Output Programs. 


a 








3. Testing of Program. 

In order to devises a logical system of tests the sim- 
plified flow chart of Fig. 2 was used. This shows the 
questions that the various sections of the subroutine 
answer. In all, there are 25 distinct paths that can be 
followed. Only four of these paths lead to 4 satisfactory 


? 


problem, i.e€., one which can be solved by simple static 


tf 


Fig. 10 lists the types of problems which were solved. The 
complete print=-outs of the solved problems are collected 
ma Appendix K. 

Problems 26 through 29 are the complements of the 
Satisfactory problems of the first 25. The "bar reactions” 
output of TESTLA was used as the "allowable bar loads" of 
TESTLB. The output of TESTLB, the "allowable load torsor*" 
Bnould be, and is, the "resultant torsor” of THSTIA. 

In problems ll and 17, where there are more than 51x 
links, the storage process is such that the last link is 
erored on top of the first. This results in a nonsense 
print-out of the structure data. In problem li, for ex- 
ample, the first link has two direction cosines of one. 

In the use of TESTLA to solve a planar problem it 
is necessary to consider the possibility of having an 
invalid problem because the load is not in the same plane 


as the structure. Since this possibility is caused by an 


ae ¥ ° . 

Torsor -~ A coined word (combination of torque and 
tensor) used to indicate a generalized resultant which 
1s composed of a couple and a force. 
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error on the part of the one setting up the problem and 
Heenou strictly a factor in the stability of the structure, 
testing for this possibility was made part of the {fnput- 
Output Program. The testing of a planar problem therefore 
requires one step beyond "STADET". The subcases of prob- 
lems 22 and 23 show the results of such tests. 

It was not possible to get an accurate estimate of 
the time reguiréd to solve a problem. It takes roughly 
one second to solve an unsatisfactory problem and about 
twice as long to solve a satisfactory one. This differ- 
ence is not in the computation time but in the time re- 
quired to write the output on the magnetic tape. In fact, 
when a series of problems is being solved, only 4 slight 
pause can be detected between the end of writing one 
answer and the beginning of the output for the next. To 
these computing and writing times must be added the time 
to compile the program. This is about 75 seconds for 


TEOTIA and about 45 seconds for TESTIB. 
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ie Conclusion. 

When this project was undertaken, it was understood 
that this subroutine would eventually be used as part of 
a much more ambitious program. The objective of providing 
a routine which decides whether the structure may be 
analyzed by the methods of simple statics has been achieved. 

For the purposes of this routine a decision was made 
to have storage for only six equivalent Links, the maximum 
number compatible with a solution obtainable from statics. 
This is the only restriction which must be removed in 
future, expanded use. The arrays KPROB (kind of problem) 
and NOSOL (no solution) were set up to be used és keys to 
select the proper methods of attack. If NOSOL‘4) is set 
to one, a redundant structure 1s involved and more power- 
ful methods must be used to effect a solution. KPROB 
indicates whether the simplification of a planar problem 
ioe valid . 

While it is probable that a bona fide planar problem 
would be set up in one of the coordinate pianes, it might 
be advisable to provide a subroutine to transform the 
coordinate system so that an oblique plane problem is 
transferred to a coordinate plane. It may or may not be 
desirable to transform the results back to the original 
coordinate system. 

The next step in the generalized analysis of structures 
by digital computers could be along one of several patns. 


Since this subroutine is concerned with only the external 
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stability of the structure, a subroutine to examine the 

internal stability needs to be written. The simplifica- 
tion of considering only space framés might yield a good 
hIrs. Approximation. 

Another path would bé to either use a code or have 
an examining type of subroutine to decide the type of 
problem and then jump to one of a set of subroutines 
designed to solve a given class of problems. 

It is felt that this is one small step in removing 


some of the drudgery from engineering désign. 


26 


a 





Bibliography 


Deresiewicz, H., Elements of Enginéering Statics, 
Columbia University Press, 1959. 


Kinney, J. Sterling, Indeterminate Structural Analysis, 
Addison-Wesley Pub. Co. Inc., 1957. 


McCracken, D. D., A Guide to Fortran Programming, John 
Wilsy and Sons, 1961 : 


Shames, Irving H., Engineering Mechanligs--Stati 
pretoece=Hail, Ines,-1959. 


Timoshenko, S. and Young, D. H., Engineering Mechanics, 
4th Edition, McGraw-Hill, 1956. 


CO OP Manual, Publication 067A, Vol. I, Control Data 
Corporation Computer Division, Dec. 1960. 


Fortran System for the Control Data 1604 Computer, 
Publication 087A, Control Data Corporation Computer 
Division, 1961. 





APPENDIX A 


Linkage Diagrams of Common Supports 


A link is a member capable of supporting only an 
axial force. It has spherical hinges at both ends so 
that it cannot transmit any moments. In the drawings 
below the link is shown thus: \ , where the dot indi- 
cates the hinge at the point of application. The other 
end of the link is to be understood to have a hinge which, 
in turn, is attached to a frame of reference. 

As discussed in Section l, Introduction, some of 
these supports are not the equivalent of a complete link. 
In these cases it is the user's responsibility to either 
modify the support or to make sure that the reaction is 


in the allowed direction. 


Smooth plane support 


Because friction is neglected, 
a linkage system which only 
opposes vertical motion is 
obtained. The link shown may 
only have a compressive reaction. 





GROUND PLANE 


Flexible cord 


GROUND The reaction is along the cord. 
PLANE Only a tensile reaction is 
allowable. 
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Ball and socket 


OCKET 


A ball and socket is the 
ae equivalent of three links 
acting at a point. 
ROUND 
PLANE 


Bui tt=-in support 







BODY BODY A built-in support is the 
vw equivalent of six links. 
ee The transparent body is 

A ea Shown supported by one of 





many stable and determinate 
GROUND arrangements of six links. 


PLANE 





APPENDIX B 


Mathematical Tests 


The examination of the system of supporting links is 
accomplished by the application of a few simple mathemati- 
cal tests to the vectors which are part of or are generated 
from the input data. 

A. The dot or inner product of two vectors A and B, 
Beauner OL which 1S zero, 1S zero if and only if A and B 


are perpendicular. 


Ee 8 = (Aj |B] cos « where a = included angle 

eee - AL BY + ALB. ae Ee lee SUDSCripis Incicete x, 
y, Z2 components of A 
and B 


Be Une ¢ross or outer product ol two vectors A and 


B, neither of which is zero, is zero if and only if A and 


Beare parallel or collinear. 
Ax Bes [A] |B] sina u where u is a unit vector per- 
pendicular to both A and 
B and with them forming 
a right-handed set 
ot a i, j and k are unit, orth- 
Eee B= | AS AL iN ogonal base vectors in the 
BBB xX, y, 2 directions, re- 
a ae”! spectively 


GC. The question of whether of not lines intersect is 
solved by examining the equations of the planes determining 
the lines. If a set of four equations détermining a pair 
of lines is inconsistent or dependent, the lines are parallel 
or intersect. The "matrix singularity test" of REAC] {see 


Appendix F) is used to find if the equations are dependent 


pLe. 





Sreinconsistent. 
If the lines intersect, the point of intersection 
can be found by solving three of the four equations simul- 


taneously. 


ful 
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APPENDIX C 
Input-Output Programs 


These control programs provide for reading in data, 
operating the required Subroutines and printing results. 
In both cases the first step is reading in the data on 
the equivalent linkages. Then "STADET" is called and the 
stability and determinateness of the support system is 
examined. The rest of the data for that problem is then 
read in regardless of the results of "STADET"”. The 
intermediate operations of the two programs will be cone 


Sidered separately. 


TESTIA 

If the structure is stable and determinate, the re- 
sultant of the forces and couples acting on the structure 
is calculated. This establishes the B vector in the matrix 
equation, Ax = B. The A matrix, that array which, when 
used to premultiply the link reaction vector, x, gives B, 
Me tound next. If there are six links, A is six by six; 
the six equations represented are independent and the values 
of the link reactions are found by subroutine REACI. 

In the case of a planar problem, "STADET" has examined 
the structure and found it satisfactory, but it remains to 
be seen whether the load is applied in the correct plane. 
Therefore, checks are made to see that the resultant force 


vector 1S in the plane of the structure and that the re- 
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sultant couple vector is perpendicular to that plane. 

If there are only three links, A is six by three; 
only three independent equations exist among the six, 
and it is necessary to solve various sets of three equa- 
tions until a consistent set is found (See Appendix F, 
REAC]). When a consistent set is found, the results are 


printed out as discussed below. 


TESTIB 

If the structure is stable and determinate, the A 
matrix is calculated. The allowable loads in the various 
links are knowns in this problem. Therefore, the x vector 
is multiplied by A giving the allowable load on the struce 
ture. In this case there is no question about the plane 


erewne resultant load as in TESIILA. 


Output 

In both programs the intermediate calculations are 
suppressed if the structure is not stable and determinate. 
The printout then consists of the input structure, a 
description of the problem type and the error. If the 
structure is satisfactory, the input for the intermediate 
calculations and their results are printed. Since the 
case of coincident support points would be an unusual 
blunder, a special error printout was not used. Since 
"STADET" is not really operated in this case, the printe- 


outs are: 
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NOT A PLANAR PROBLEM and 


STRUCTURE UNSTABLE LINES OF ACTION OF SUPPORTS 
INTERSECT A LINE 


(See problem 1.) 
Both programs reset conditions so as to solve the 
next problem and then continue until all solutions have 


been made. (See Appendix H for problem format. ) 
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APPENDIX D 


SUDrOUtI nee.) . 


The subroutine evaluates the cross product of two 
vectors expressed in rectangular coordinates. A Fortran 
expression is written for each term of the expansion of 


CROSS. 
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APPENDIX E 


Subroutine RES1 


The subroutine finds the resultant of a system of 
forces and/or couples by first adding up the rectangular 
components of either and then calculating the effect, if 


any, of the forces on the resultant couple. 
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LoOnce 
ee ry and r, are the rectan-~ 
gular components of r 

ae Pa = es F is the vector force 
Ree By and Eo are the réctan- 
gular components of F 
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APPENDIX F 
Subroutine REACL 


This subroutine is a modification of the GAUSS2 sub- 
routine used by C. B. Bailey“ in his program for solving 
Simultaneous linear equations. Whereas Bailéy was solving 
the matrix equation, Ax = B, with a maximum of 50 vectors 
B for each coefficient matrix A, nowhere in this program 
is more than one B vector involved with a given A. 

The routine locates the largest first column coef- 
ficient and, if necessary, exchanges the row containing 
this coefficient with the first so the largest element is 
Aji: Multiples of row one are subtracted from the other 
rows to make the column one coefficients below the first 
zero. Next the column two coefficients below the first 
are examined to find the largest and rows are exchanged, 
if necessary, as before. This Gaussian elimination process 
is continued until all the elements below the main diagonal 
are zero. 

At each reduction step, the value of the diagonal 
element is compared with zero. A very small diagonal 
element indicates ill-conditioned equations and causes 
an error output. This is the "matrix singularity" test 
used in CONCURRENT1. If the matrix is non-singular, a 
back- solving process is used to find the unknowns. 


“co. B. Bailey, F2 UTEX LINEQN, CO OP Manual, Dec., 1960. 
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APPENDIX G 


Subroutine PLANE 


Given a vector with direction cosines 1, m and n 
and one point, the problem is to find a pair of planes 
whose intersection is parallel to the vector and passes 
through the point. 

First, the subroutine checks to see if the vector 
is "nearly" parailel to a coordinate axis. "Nearly" is 
defined as having a direction cosine greater than 0.8, 
Z.€., lying within about 37° of the coordinate directions. 
Bi the vector is "nearly" parallel to one axis, it is 
crossed with vectors along the perpendicular axes. 

If the vector is not "nearly" parallel, a check is 
made to see whether it lies within about 6° of a coordinate 
plane, i1.¢€., the direction cosine is less than O.1. If 
it is close to a coordinate plane, one of the crossing 
vectors is chosen perpendicular to that plane. 

If neither of the above tests is met, the given 
vector is crossed with vectors in the x and z directions. 

There are now two cross product vectors of the form: 


X, = AL + BU + Ck 


Since these vectors are perpendicular to the pianes con- 
taining the vectors which formed the cross products, the 
dot product of X, and an arbitrary vector in one of the 
planes is zero. An arbitrary vector is of the form: 


(x=X, jae a (y-y, yal a (2-2,)k 


60 





where (x,,y,,2,) is the known point. The dot product has 
the form: 
Mixa) ota yay) et 2524) = 0 
This may be rearranged to: 
Oe ar a eG 2. ae 
moere ; ID oy WheKee? oye. a. 80 Ze 
The coefficients A, B, C, and D are stored in the array 


PLA. 
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APPENDIX H 


Problem Format 


TESTLA 

Fortran Format Symbol(s) 
5 NPROB 
Hollerith 

ICS NBAR 
BrlO.6, 3 Ply. 7 Noes 

fee 12 NFO , NCO 
6F12.6 FORCE, CORD 
Lifes 7 COUPL 


Input data 
Total number of problems 
Problem title - up to 80 
characters 


Number of links 


Direction cosines, coordi-~ 
nates of links 
Number of forces, couples 


Force components, points of 
application 


Couple components 


Repeats from the Hollerith title for each problem. 
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TEST1B 


Fortran Format 
J 
Hollerith 


5 
plo.8, 3F14.7 


Orie .6 


Symbol(s) 
NPROB 


NBAR 
VB,B 


BFORCE 


Input Data 

Total number of problems 
Problem title = up to 80 
characters 

Number of links 

Direction cosines, coordi- 
nates of links 


Allowable link loads 


Repeats from the Hollerith title for each problem 
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KPROB = 


1 


a 
3 
4 


NOSOL = 


APPENDIX I 
STADET Flags 


Kind of Problem 

Support points are on a line 
Support points are in a plane 
Reaction matrix is singular 
Planar problem (if 4 is 1) 

Not a planar problem (if 4 is 0} 


Resultant load is not in plane of structure 


No solution because: 

Structure unstable - lines of action of supports 
are concurrent 

Structure unstable - lines of action of supports 
are parallel 

Structure unstable - insufficient number of 
supports 

Structure indeterminate - excess number of supports 
Structure unstable - lines of action of supports 
intersect a line 

Structure unstable - supports are in parallel 


planes 
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APPENDIX K 


Problem Types 


In this appendix "bars" and “supports” are used as 


synonyms for "links". See Section 3, Testing of Program, 


for a discussion of the problems. 
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APPENDIX L 


Detimncdion of “Zero” 


The computer must convert any non-integer number to 
an approximate binary form which will fit its word length. 
The conversion and length limitation allow an accuracy of 
about ten decimal digits. With each calculation round-off 
errors accumulate. Whereas the flow charts ask if a 
quantity is zero, what must really be considered is 
whether the quantity is less than the maximum expected 
round-off error. If there are no conflicting requirements, 
the maximum round-off error may be estimated and, with 
perhaps some leeway, defined as "zero". However, where 
different orders of magnitude may be involved in the same 
problem or in subsequent problems or where it is advisable 
to have a large tolerance for one calculation and a small 
tolerance for another, the situation becomes far more 
complicated. One cannot define the "zeros" in terms of 
the number of significant figures in an expression; 4 
specific magnitude is required. 

An attempt was made to tie the size of the "zero" to 
the size of the structure. It was felt that barring un- 
reasonable complications, such as having a very small 
structure at a great distance from the origin, the system 
of averaging the coordinates of the points of application 
and multiplying this average by the computer accuracy 


would give a suitable zero. However, instead of being 


ii aeae 


os 





able to use ee times the average coordinate, it was 
necessary to use a "zero" 10,000 times as big to get 
proper answers. Further investigation revealed that 
the difficulty probably involved the use of subroutine 
REACL in CONCURRENTI1. 

In subroutine REAC] a small value of "zero" allows 
small diagonal elements which, during the back-solving 
process, generate large coordinates of the point of 
intersection. When these large values are then used 
in CONCURRENTL, the round-off error becomes much greater 
than the "zero" defined above. Unsuccessful efforts were 
made to modify the "zero" obtained above in a logical way 
for use in CONCURRENT1 and REAC1. 

6 


A multiplier of 10 ~ gives proper answers, at least 
for this set of problems, but there is no other justifi- 
cation for its use. The author, regretfully, did not 


have sufficient time to pursue this investigation. 
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